Space and Scale Localization in a Model of a Direct Transition 

to Spatiotemporal Chaos 



Guoming Xiong* 
Department of Physics, Lehigh University, Bethlehem, PA 18015 

J. D. Gunton 

Department of Physics, Lehigh University, Bethlehem, PA 18015 

Haowen Xi 
Department of Physics and Astronomy, 
Bowling Green State University, Bowling Green OH 43403 

Raul Toral 

Institute Mediterraneo de Estudios Avanzados (IMEDEA), CSIC-UIB 
Campus Universitat de les Illes Balears, E-07071 Palma de Mallorca, Spain 

February 3, 2008 

Abstract 

We present the results of a wavelet-based approach to the study of the chaotic dynamics of a 
one dimensional model that shows a direct transition to spatiotemporal chaos. We find that the 
dynamics of this model in the spatiotemporally chaotic regime may be understood in terms of 
localized dynamics in both space and scale (wavenumber). A projection onto a Daubechies basis 
yields a good separation of scales, as shown by an examination of the contribution of different 
wavelet levels to the power spectrum. At most scales, including the most energetic ones, we 
find essentially Gaussian dynamics. We also show that removal of certain wavelet modes can be 
made without altering the dynamics of the system as described by the Lyapunov spectrum. 
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1 Introduction 



A major goal in the study of temporal chaos in spatially extended systems |l], ||, |j (spatiotemporal 
chaos or STC) is to find a statistical description of the behavior of a particular dynamical system in 
the limit of large length and long time scales. This is analogous to a hydrodynamic description of a 
system of microscopic particles satisfying classical mechanics. Some progress has been made in this 
regard in recent years, including constructing a long wavelength, long time theory of the Kuramoto- 
Sivashinsky (KS) equation . In this reference, the authors obtained an effective stochastic equation 
which belongs to the Kardar-Parisi-Zhang (KPZ) university class in the hydrodynamic limit. This 
was obtained by incorporating the chaotic dynamics of the small KS system in a coarse-graining 
procedure. The basic premise of their approach is that the spatiotemporal chaos of a large system can 
be understood in terms of the chaos observed in mutually coupled, small systems. Two other recent 
studies |(| of the Kuramoto-Sivashinsky equation used a wavelet decomposition to characterize 
its spatiotemporal chaos. The first investigation led to results similar to those of reference 0], 
suggesting a statistical description of a group of identical short length subsystems, slowly driven via 
interactions with the larger scales. More precisely, in || it was shown that an effective equation 
could be obtained and consistently approximated by a forced Burgers equation, for scales far from 
the cutoff between small and large wavelengths. This work was extended in || where the authors 
found that projecting onto a spline wavelet basis enabled a good separation of length scales, with 
each having its own characteristic dynamics. At large scales they found essentially slow Gaussian 
dynamics, which can be understood in terms of local events. The results are also consistent with 
the picture of weakly interacting small subsystems and so-called "extensive chaos" (in which the 
Lyapunov dimension is proportional to L d , where L is the linear system size and d the dimension 
of space) . The authors also discussed various correlation lengths and demonstrated the existence of 
a spatial interaction length, which provides a limit on how much one may limit spatial interactions 
without changing the dynamics significantly (and hence limits how small a system one can use). 

Motivated by the success of the wavelet decomposition of the KS model, we have undertaken 
such a decomposition for another one dimensional model, the so-called Nikolaevskii equation |?], ||, ||] . 
This model was originally proposed to describe the propagation of longitudinal seismic waves. It 
was subsequently shown to exhibit a direct transition from a spatially uniform, stationary state to 
a spatiotemporally chaotic state as a control parameter was varied Q . The spatiotemporal chaotic 
behavior of this system (which has been named as soft mode turbulence or SMT) is different from 
that of the KS equation as a result of the existence of an additional continuous symmetry in the 
model (beyond the conventional symmetries of space and time translation invariance). This model 
is particularly interesting in that there are two control parameters, e and the system size L, such 
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that one can study the transition to chaos as a function of e in terms of power law behavior, scaling, 
etc. In this sense it is a richer model than the KS equation and many results have been obtained for 
it, including a calculation of the Lyapunov exponents, Lyapunov dimension and Kolmogorov- Sinai 
entropy for several values of e in the limit of large L such that extensive chaos holds jlO| . It has also 
been shown that the distribution function for the order parameter is Gaussian for large wavelengths 
and large times. In this paper we carry out a wavelet decomposition similar to that in reference 
||, using a Daubechies basis. We find that the most energetic modes have a Gaussian distribution 
for this choice of basis. We also calculate the Lyapunov spectrum and find that one can remove a 
certain set of the modes without altering this spectrum, suggesting that in this sense one can obtain 
a more minimal description of spatiotemporal chaos for this model. 

The outline of the paper is as follows: in section II we define the model and the wavelet decompo- 
sition scheme. In section III we present the main results of our analysis. This includes a calculation 
of the energy distribution (power spectrum) and the probability density functions for the wavelet 
coefficients. It also includes a calculation of the effects on the Lyapunov spectrum resulting from 
removing various modes from the dynamics. In section IV we present briefly the conclusions of this 
work. 



2 Model and Wavelet Decomposition 



The Nikolaevskii model is defined by the partial differential equation 

9v d 2 , d 2 , 2l dv . 

with < x < L and periodic boundary conditions, where e and L are two control parameters for 
the model. In the present paper we study the model for e = 0.5 and L w 158; these values are large 
enough to see well developed STCjTo|. 

The general form of a discrete wavelet decomposition for a field v(x,t) can be written as 

"fe^EE^w^w ( 2 ) 

j k 

where the set tyj k (x) = 2 5 ' 2 ^(2 5 'a; — k) form an orthogonal basis in the sense that 



(tfjfc | 4v fe <) = / * jk (x)V fk '{x)dx = S j}f 5 k , k , (3) 

J — OO 

The indices j, k are integers which we specify below. The function *ff(x) is called the wavelet function 
(or "mother" function) and the wavelet coefficients can be obtained as 

a jk (t) = (v | V jk ) = / v(x,t)$ jk (x)dx. (4) 
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One obtains equations of motion for the wavelet coefficients by substituting Eq.(g) into Eq.(Q). 
Those equations of motion can then be numerically solved in order to obtain the time evolution of the 
wavelet coefficients and, hence, information about the chaotic behavior of the model. Alternatively, 
one can solve directly Eq. (jl|) and use Eq.(Q) to obtain the wavelet coefficients. This latter approach 
is the one we follow in this paper: at the required times we decompose the resulting solution for 
v(x,t) in terms of the Daub4 orthogonal basis set. The Daub4 is the simplest of a wavelet family 
named DaubA" constructed by Daubcchies Jll[ where A ranges from 4 to 20. Based on this wavelet 
family a very effective algorithm has been developed pl| . Thus by this method we obtain the time 
dependence of the wavelet coefficients a,jk(t). 

The Nikolaevskii equation ([!]) is integrated numerically using a version of the pseudo-spectral 
method combined with a fourth-order predictor-corrector integrator. (The details of this method 
will be given elsewhere). We use a mesh size 5x — 0.31 and a time step 5t = 0.01 in the numerical 
scheme. The length L of the system is L — N5x, where N is the total number of points of the 
system. It is convenient to choose N = 2 J+1 (J is an integer which represents the largest level 
of the wavelet) in order to apply Daub4 in an efficient pyramidal scheme. In our case we usually 
choose J = 8soA = 512 and L = 158.72. The total integration time in our simulation goes from 
t = 600,000 to t = 1,100,000 (in the dimensionless units of equation ([j])), depending on which 
quantity we calculate. Usually, once we have reached the chaotic state, we perform time averages of 
the quantities of interest. 



The wavelet decomposition using Daub4 is carried out in a pyramidal scheme such that at level 
j there is a total of 2- 7 coefficients: ajk{t) with k = 1, . . . , 2 J and j = 0, . . . , J, except for the coarsest 
level j — where there are 2 coefficients, namely a 00 (t) and aoi(i), instead of just one. To be more 
precise in the notation, we define the wavelet contribution at level j as 

Vj(x,t) = y^^ajk(t)^jk(x), (5) 

k 

such that the field is a sum of all its wavelet contributions, v(x, t) = J^j v j( x i Notice that as the 
wavelet level becomes coarser, i.e reducing j by one, the number of wavelet coefficients contributing 
to Vj(x, t) is therefore reduced by a factor of two and that the total number of wavelet coefficients 
ajk{t) equals the number of points in the system, N = 2 J+1 . 



3 Wavelet Analysis: Results 
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3.1 Energy Distribution 



In the language of signal processing, the total energy of a signal field is a conserved quantity. For 
our case, the energy for the chaotic field v(x,t) can be defined as 

E(t) =L-\v | v) = L~ l [ v 2 {x,t)dx (6) 

Jo 

and it can be written as the sum of the energies at each wavelet level, E(t) = X^ e .?'W with 
ej(t) = L^lyVj | vj). Fig. (la) shows the energy distribution per wavelet level for the Nikolaevskii 
model, averaged over time. One sees that the energy of the field is mainly concentrated at wavelet 
level j = 5. As a comparison, we also show the structure factor S(q), which can be thought of 
as the energy distribution in Fourier space, in Fig. (lb). It is easy to show for this model || 
that the unstable modes (within linear stability analysis) in Fourier space are located between 
<7i.2 = (1 ± -v/e) 1 / 2 , with the most unstable mode at q m — 1.0, consistent with the results shown in 
Fig. 1(b). In our case with e = 0.5 it is qi ~ 0.54 and qi ~ 1.31, so with L = 158.72, the most 
energetic mode should be at n m = q m L/2ir ~ 25.26 and the smallest and largest unstable modes are 
at ni = qiL/2n ~ 13.67 and n-i = <hl^ — 33.00 respectively. Therefore there should be about 20 
complex or 40 real unstable or marginal modes concentrated in the neighborhood of the 50th real 
mode. Notice that wavelet levels — 4 contain a total of 32 modes (corresponding to 32 wavelet 
coefficients). It is easy to understand, therefore, that the energy peak is located at wavelet level 
j = 5. On the other hand, if one compares Figures (la) and (lb), we can also see that for small 
q there is no similarity between the two energy distributions, implying that for small q the overlap 
between different wavelet levels is relatively strong. 



3.2 Temporal Behavior and Probability Distribution of Wavelet Coeffi- 
cients 

Figure 2 shows the temporal behavior of some wavelet coefficnts a,jk(t). At each level j we show the 
intermediate wavelet coefficient with fc = 2 J . Notice that, from the very early times, there is a 
chaotic behavior in the temporal evolution. In Fig. 3 we show the probability density 

functions (PDFs) for the wavelet coefficients from level j = 1 to j = 7, averaged over time. 
Namely, for fixed j we sum over all k for the a,jk coefficients to obtain a mean value; we repeat this 
over all times sampled and plot these mean values to obtain the PDF. The PDFs for j = and j = 8 
are not shown here, since they are just Dirac-delta functions with spikes at zero. The delta function- 
like behavior at the largest wavelet level j — S simply means that there is a strong dissipation at 
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large wavenumbers in the Fourier spectrum (cf Wittenberg and Holmes ||). (In Fig. 3 we can 
already see a similar behavior developing at the level j — 7.) The spike at j = 0, however, only 
implies that the coarsest part is not suitable for describing the field and its two components compete 
with each other everywhere to yield a zero average value for the amplitude. We have tried to fit 
these PDFs to a Gaussian distribution. In Fig. 3 we see that for j — 1, 7 the PDFS are essentially 
Gaussian. (We do not think the small deviations from Gaussian are statistically significant.) 

The wavelet decomposition of the STC in this model allows us to examine the chaotic behavior 
of the system at different spatial scales. In general, as can be seen somewhat in Fig. la), the large 
j levels correspond to large wavenumbers or small spatial scales and the small j levels correspond 
to small wavenumbers or large spatial scales. In this section we try to confirm this in more detail 
by showing the structure factors which result from removing various wavelet levels. 

In Fig.(4a-f) we show the structure factor for some of these "reduced level" systems, with the 
full structure factor for the Nikolaevskii model shown for comparison. Fig. (4a) shows the result 
of removing wavelet level j = 8, for which the new structure factor is almost identical with the 
exact one. This shows that removing the smallest spatial scale in the wavelets (or, equivalently, 
the shortest wavelengths) leaves the system invariant. Fig. (4b) shows the result of removing levels 
j = 7 and j = 8. In this case, there is a small difference between the two structure factors for 
large wavenumbers, indicating that the large j levels of the wavelets only contribute to the short- 
wavelength dynamics (i.e. the fast dissipation) of the model. Fig. (4c) shows the result of removing 
levels j = 6, 7 and 8. In this case there is a significant difference between the two structure factors 
near the peak. This is not surprising since j = 6 is the second most energetic level in the model. Fig. 
(4d) shows the results of removing levels j = 0, 1 and 2. The resulting structure factors shows that 
these levels are responsible for the small peak in the structure factor near wavenumber q w 0. Fig. 
(4e) shows the role of the most energetic levels in the model, i.e. j = 4, 5 and 6. We see that these 
levels are responsible for the dominant part of the structure factor. Furthermore, we also see that 
without level j = 3, the small peak near q ss becomes even smaller in comparison with Fig.(4d). 
Finally, Fig. (4f) shows the result of removing all levels with j > 2. As one would expect, one only 
has the peak near small q. If one considers all of the above results, one obtains a clear picture of 
the scale localization involved in the wavelet decomposition. 

We end this section with a discussion of the relationship between wavelet decomposition and 
the renormalization group transformation. It is clear that removing the large- j wavelet levels is a 
procedure that essentially removes the short-wavelength part of the structure factor in Fourier space. 
In this sense it is like a renormalization group transformation in which one integrates out the short 
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wavelength degrees of freedom to obtain a renormalized description of a system. The particular 
Daub4 algorithm we use here is a kind of decimation transformation in that at every step removing 
a wavelet of the largest remaining level means that one has only half as many spatial points as before 
the removal to describe the (chaotic) field. The other set of values of the field is "filtered out" , with 
the value of the field at a given new (remaining) space point being a kind of average of the field at 
four of the neighboring points of the original set. Other types of wavelet decompositions correspond 
to different kinds of transformations. 



3.3 Lyapunov Spectrum 



We conclude by examining the effect of removing wavelet levels on the Lyapunov spectrum of the 
system. This provides us with the most detailed understanding of the contribution of various levels 
to the chaotic dynamics of the system. Fig. (5a) shows the Lyapunov spectrum obtained after 
removing the j = 8 level, while Fig. (5b) shows the Lyapunov spectrum obtained after removing 
j = 7 and j = 8. The Lyapunov spectrum for the Nikolaevskii model is also shown for comparison. 
In Table 1 we show the Lyapunov dimension D and the Kolmogorov-Sinai entropy H for these 
different cases. It seems that to a very good approximation the chaotic behavior of the system 
does not depend on the j = 8 level. Thus a simplified statistical description of the Nikolaevskii 
equation can be obtained by excluding the j = 8 level, without altering the chaotic dynamics for 
the system. However, it is also clear from inspection of Table 1 that removing any additional levels 
significantly alters the behavior. Thus all these levels play an important role in the spatiotemporal 
chaotic dynamics of the system. 



4 Conclusion 



We have shown that the dynamics of the one dimensional Nikolaevskii equation in the spatiotem- 
porally chaotic regime may be understood in terms of localized dynamics in both space and scale 
(wavenumber). Specifically, a projection onto a particular Daubechies basis (Daub4) yields a good 
separation of scales, as shown by an examination of the contribution of different wavelet levels to 
the power spectrum. At most scales, including the most energetic ones, we find essentially Gaussian 
dynamics. Perhaps most importantly, we also found that removal of certain wavelet modes can be 
made without altering the chaotic dynamics of the system as described by the Lyapunov spectrum. 
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Many different length scales have been proposed for the description of spatiotemporal chaos g]. 
These include the usual correlation length £2 for the (two point) order parameter and the dimension 
correlation length £5, obtained from the Lyapunov dimension D(L) as £5 = limL^^U 1 / D{L) . Here d 
denotes the space dimensionality and L the linear dimension of the system. We have studied both of 
these lengths as a function of the control parameter e for this model. The results of these calculations, 
as well as for other quantities characterizing the spatiotemporal chaos will be presented elsewhere 
|P^| . Since it is difficult to calculate the Lyapunov dimension for high-dimensional systems, Zoldi 
and Greenside have proposed using a Karhunen-Loeve dimension D^ld, defined by the number 
of eigenmodes in a proper orthogonal decomposition necessary to capture a given fraction / of 
the total energy. From this dimension one can define a Karhunen-Loeve correlation length ^kld- 
However, for a translationally invariant system, such as the Nikolaevskii equation with periodic 
boundary conditions, the Karhunen-Loeve eigenmodes are Fourier modes. Thus the Karhunen- 
Loeve dimension length for any / can be computed directly from the power spectrum S(q) and thus 
contains no more dynamical information than £2 ■ As Wittenberg and Holmes have pointed out || all 
these lengths (with the possible exception of the dimension correlation length £g) are measures only 
of spatial disorder and thus yield no information about the temporally chaotic dynamics responsible 
for the disorder. As a consequence, Wittenberg and Holmes introduced another length scale, a so- 
called dynamical interaction length. This is a length scale l c such that if one deletes interactions for 
length scales greater than l c in the wavelet Galerkin projection of the model equation of interest, 
one alters the chaotic dynamics of the model (i.e. changes the Lyapunov spectrum). Although this 
is in principle a very interesting length scale to study, it involves a numerical calculation which for 
our model is computationally expensive and beyond our current resources. As a consequence, we 
leave this important issue for future investigation. 

Finally, as mentioned in the introduction, Wittenberg and Holmes [[| performed a wavelet de- 
composition of the Kuramoto-Sivashinsky model, using an orthogonal spline basis. They found that 
the probability distribution functions for the most energetic modes were non-Gaussian. This differs 
from our results for the Nikolaevskii model, in which the distribution functions for the most ener- 
getic modes are Gaussian (to a very good approximation and excluding the peak near zero). To see 
whether this difference is significant or simply results from a different choice of basis functions in 
the two studies, we have calculated the PDFs for the KS model using the Daub4 basis. We found 
that with this choice of basis the PDFs for the most energetic modes are in fact Gaussian. We 
find non-Gaussian behavior only at the small- j levels. Thus it would seem that any non-Gaussian 
behavior of the wavelet coefficients depends on the basis chosen for the wavelet decomposition. 
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Wavelet levels 


Lyapunov dimension 


KS entropy 


including j=0-8 


52.31 


5.27 


including j=0-7 


51.89 


5.15 


including j=0-6 


51.40 
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including j=0-5 
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including j=0-4 
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including j=4-8 


57.10 
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including j=4-6 


57.00 


7.95 


including j=0-2,4-8 


55.13 
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Figure 1: (a.)The energy distribution per watelet the field v(x, t) satistying the Nikolaevsky equation 
(1) with e = 0.5 in a finite length L ps 158. The results shown come from a numerical integration 
and have been averaged over time. Notice that the peak is located around the 5th wavelet level. 
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Figure 2: The temporal behavior of the wavelet coefficients ajk for different levels j and a selected 



intermediate value of k. 
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Figure 3: The probability distribution function for the wavelet coefficients ajk for levels j = 1 to 
j = 7, averaged over time. 






Figure 4: The structure factor which resilM from removing different wavelet levels (dashed line) 
compared with the full structure factor (solid line). The levels in each figure are: (a) level j = 8 
has been removed, (b) levels j = 7,8 have been removed, (c) levels j = 6,7,8 have been removed, 
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Figure 5: The Lyapunov spectrum obtained including all the wavelet levels (solid line) compared 
with that obtained removing level j = 8 (top figure) and levels j = 7,8 (bottom figure). 



